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,—1 Abstract 

CIh It is well known that i?-symmetric models dramatically alleviate the SUSY flavor and 

o, CP problems. We study particular modifications of existing i?-symmetric models which 

QJ share the solution to the above problems, and have interesting consequences for electroweak 

baryogenesis and the Dark Matter (DM) content of the universe. In particular, we find 
that it is naturally possible to have a strongly first-order electroweak phase transition 
while simultaneously relaxing the tension with EDM experiments. The i?-symmetry (and 
Q^^ its small breaking) implies that the gauginos (and the neutralino LSP) are pseudo-Dirac 

fermions, which is relevant for both baryogenesis and DM. The singlet superpartner of 
the U{1)y pseudo-Dirac gaugino plays a prominent role in making the electroweak phase 

J I transition strongly first-order. The pseudo-Dirac nature of the LSP allows it to behave 

(^ similarly to a Dirac particle during freeze-out, but like a Majorana particle for annihilation 

^~^ today and in scattering against nuclei, thus being consistent with current constraints. 

. . Assuming a standard cosmology, it is possible to simultaneously have a strongly first-order 

_^ phase transition conducive to baryogenesis and have the LSP provide the full DM relic 

S^ abundance, in part of the allowed parameter space. However, other possibilities for DM 

H also exist, which are discussed. It is expected that upcoming direct DM searches as well as 

neutrino signals from DM annihilation in the Sun will be sensitive to this class of models. 

Interesting collider and Gravity-wave signals are also briefly discussed. 
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1 Introduction 

Although the Standard-Model is an excellent description of Nature up to energies of around a 
hundred GeV, it fails to explain two of the most important mysteries of our Universe - the nature 
of Dark Matter (DM) and the origin of the matter-antimatter asymmetry. F] Many candidates 
for Dark Matter - axions, gravitinos, WIMPs, etc. exist; however the most popular among 
them is the WIMP, which can provide the observed relic abundance from a thermal freeze-out 
mechanism in a large region of parameter space. Moreover, such a particle also arises within 

^We have nothing to say in this paper about the biggest mystery - the tiny value of the Cosmological Constant. 



many schemes of beyond-the-Standard-Model physics such as low-scale super symmetry, extra 
dimensions, etc. 

Similarly, many popular mechanisms for generating the matter- antimatter asymmetry exist. 
These can be broadly divided into those which depend on physics at very high scales and those 
which depend on physics at the electroweak (EW) scale and below. Models in the first class 
include models of GUT baryogenesis [1] and leptogenesis |2]. On the other hand models in 
the second class include those of electroweak baryogenesis [3], TeV-scale leptogenesis [1] and 
Affleck-Dine baryogenesis [S] . There also exist models which propose to explain both the baryon 
asymmetry of the universe (BAU) and the DM abundance as an asymmetry in a conserved 
charge [6]- [11], [12], although they typically involve a very weakly coupled "hidden" sector. 

Electroweak baryogenesis is attractive since it only depends on physics in the visible sector at 
the electroweak scale; hence the mechanism is completely testable (at least in principle) at col- 
liders probing energies of order the electroweak scale, or astrophysical/cosmological observations 
probing temperatures (or times) of order the electroweak scale. Supersymmetric extensions of the 
Standard Model provide an elegant solution to the hierarchy problem as well as a Dark Matter 
candidate in the form of the lightest supersymmetric particle or LSP (assuming i?-parity con- 
servation). However, within the minimal supersymmetric Standard Model (MSSM), the allowed 
region for a strongly first-order phase transition, which is one of the important requirements for 
electroweak baryogenesis, is severely constrained [13] and will soon be tested by experiments 
(see [H] for a recent review). In addition, there is tension between having enough CP violation 
to produce the required baryon asymmetry and constraints from electric dipole-moment (EDM) 
searches [15]. It is, therefore, worthwhile to study extensions of the MSSM with additional terms 
in the tree-level scalar potential which could help in making the transition strongly first-orderjj 
The viability of singlet extensions of the MSSM with regard to electroweak baryogenesis has 
been widely studied [18] [19]. In particular, Ref. [19] studied a singlet extension of the MSSM 
called the nMSSM, in which it was shown that it was possible to have a strongly first-order 
phase transition as well as generate the correct relic abundance for the lightest neutralino. The 
lightest neutralino in that case is, however, not completely stable but has a lifetime longer than 
the age of the Universe. 

In this work, we explore a different extension of the MSSM which can give rise to electroweak 
baryogenesis and a WIMP DM candidate. Assuming standard cosmology, the freeze-out abun- 
dance of the WIMP can account for the full DM abundance in part of the allowed parameter 
space. It is also possible for the WIMP to form an 0{1) fraction of the DM, or a negligible frac- 



^ Alternatively, it has been argued that new fermions, strongly coupled to the SM, can also strengthen the 
EWPT, and lead to a EWBG/DM connection [T^]. Ideas for generating the baryon asymmetry at the TeV scale 
without a strongly first-order EWPT have also been proposed [T7] . 



tion of DM, in other regions of the allowed parameter space. The WIMP DM candidate in the 
framework is completely stable. In addition, the tension with EDM searches is largely relaxed or 
inexistent. A crucial feature of the model is the existence of an (approximate) ?7(l)/j- symmetry, 
leading to the following important consequences: 

• The i?-symmetry leads to new interactions in the scalar potential which help make the 
phase transition strongly first-order. 

• Since global symmetries are expected to be violated by gravitational effects, "small" R- 
breaking effects are included. The i?-symmetry can thus be thought of as an accidental 
symmetry of the low-energy theory. As we will see, the small i?-breaking helps in achieving 
many phenomenologically desirable features at the same time: 

a) Viable down-type fermion masses in some realizations (others do not require R-breaking). 

b) Evade stringent direct- detection bounds from XENONIOO. 

c) Obtain a relic abundance consistent with observations in part of parameter space. 

• The CP-violation required for baryogenesis may have two origins: a) There are explicit CP 
phases in the parameters of the scalar potential, resulting in complex, spacetime-dependent 
vevs, inducing CP-violating phases in the mass matrices for the charginos and neutralinos, 
or b) All parameters in the scalar potential are real and CP violation arises from the phases 
in the i?-symmetry breaking Majorana masses of the charginos and neutralinos. In both 
cases, the approximate i?-symmetry suppresses the contribution to EDMs relative to that 
in the non i?-symmetric case. 

The plan of the paper is as follows. In Section |2] we present the relevant features of the 
i?-symmetric scenarios, and discuss a number of realizations and theoretical considerations. In 
Section [3| we discuss the zero-temperature properties of such models, while in Section |4] we 
analyze the finite-temperature potential and explain the origin of the strongly first-order phase 
transition. In Section Isl we discuss several aspects related to DM (relic density, direct and 
indirect searches), and its connection to the first-order electroweak phase transition (EWPT). In 
SectionlTlwe briefly comment on the generation of the baryon-asymmetry of the Universe (BAU), 
and discuss the reasons that allow to satisfy EDM constraints, even with sfermion masses at 
around a TeV. Finally, we consider other signatures in Section [8} and conclude in Section [9j We 
also include two appendices with technical details. 

2 Features of the Model 

In this section, we describe the structure of the model at the electroweak scale in detail. An 
i?-symmetry is well motivated from a theoretical point of view. In the global limit, having an 



exact i?-symmetry allows one to construct simple models of dynamical supersymmetry breaking 
in accord with the Nelson-Seiberg theorem [20] ■ A unique anomaly-free i?-symmetry also arises 
in a superconformal field theory at the superconformal fixed point, which can be found by a- 
maximization [21]. One of the first phenomenological models with i?-symmetry was written 
down in [22] . 

Having an i?-symmetry gives rise to many interesting phenomenological features. For ex- 
ample, Majorana gaugino masses, trilinear A parameters, and the -B/i parameter are forbidden. 
However, gaugino masses of the Dirac type are allowedjj As a result, this leads to a significant 
suppression of flavor and CP-violating effects relative to the MSSM for 0{1) flavor-violating soft 
scalar masses and phases [24j. Although the minimal i?-symmetric spectrum does not give rise 
to gauge coupling unification, many scenarios for adding additional matter have been proposed 
which could help unify the couplings |23l |25l |26l [27] . 

The minimal i?-symmetric spectrum is shown in Table 1 with the possible i?-charge assign- 
ments. The origin of these assignments can be understood as follows. The i?-charges of the 
adjoints {S,T,0} are chosen to ensure that a Dirac gaugino mass term can be written down 
with a D SUSY-breaking spurion. The simplest i?-charge assignment of the Q, t/*^, and H^ super- 
fields, so that the up-type Yukawa couplings are allowed, is given bjy-R[(5, U'^] = l,R[Hy\ = 0. 
This leaves the i?-charges of the D^, -E^, L and Hd superfields. A number of approaches can be 
followed about the -R-charges of these fields and the associated couplings: 

• I. One approach to down-type fermion masses was taken in ^27i| with R[H(i\ = 2, R[D'^, E'^] = 
—1, and R[L] = 1. This allows down-type Yukawa couplings in the superpotential, at the 
expense of making the i?-symmetry anomalous with respect to SU{3)c- However, since the 
R symmetry is best thought of as an accidental symmetry, this does not appear to be an 
issue. There are two ways of generating viable down-type fermion masses in this scheme. 
Small i?-breaking effects could arise from Majorana gaugino masses and/or A parameters, 
which could generate down- type fermion masses at 1-loop [20] ■ Alternatively, /^-breaking 
could arise from a Bfi term, leading to a small vev for H^ and providing viable down-type 
fermion masses via the usual holomorphic Yukawa couplings. 

• II. A model with R[L, D^,E^] = 1 and R[Hd\ = 2 is considered in [28], which makes the 
i?-symmetry anomaly free. In this case, down-type Yukawa couplings are forbidden in the 
superpotential, but down-type fermion masses can still be generated by supersymmetry- 
breaking terms in the Kahler potential. Generating viable fermion masses and an elec- 
troweak scale /z-term then implies that the supersymmetry breaking scale is quite low. 



•^Dirac gaugino masses can also be motivated from "supersoft" supersymmetry breaking in which the gauge 
sector has N =2 supersymmetry [23] . 

^These terms must be ahowed without suppression since the top mass arises from these. 
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Table 1: Quantum numbers of the superfields in the modeL The four entries for the i?-charges of the superfields 
D"^, E'^, Hd correspond to their values in approaches I-IV described in the main text. The fields {Ru,Rd} and 
{H[^,H'^i} only exist in approaches III and IV, respectively. 



Also, [28j considered a version of the model without the singlet S. 

• III. Both up-type and down-type fermion masses are allowed if R[Hd] = 0, and R[L, D^, E^] = 
1. However, terms like ^H^Hd and XgSHuHct which generate Higgsino masses are forbidden. 
Viable Higgsino masses can however be generated by having additional vector-like doublets 
Ru and Rd with i?-charges R{Rui Rd\ = 2, which also makes the i?-symmetry anomaly free 
[211 [30]. 

How does our approach compare with those above? The two main issues of interest in this 
paper are the existence of a strongly first-order electroweak phase transition (EWPT), and of 
a WIMP DM candidate which could provide all or some of the DM relic abundance and have 
a signal in DM experiments. It turns out that within the i?-symmetric setup, the former can 
be naturally achieved in the presence of the superpotential term XgSHuHd and the potential 
term tsS + h.c. (as discussed in Section [3]). This implies that R[Hd] = 2. The requirement of 
a strongly first-order EWPT, therefore, seems to suggest that we either follow approach I, or 
approach II (with the singlet field taken into account). However, even with the singlet, approach 
II requires very low-scale supersymmetry breaking which typically gives rise to an extremely 
light non-WIMP gravitino DM candidate -not desirable from our point of view. Approach HI 
does not work within our framework since it assumes R[Hd\ = 0. 

However, it is possible to use an approach which is a combination of that in II and HI, and 



is consistent with our requirements. We call it approach IV and briefly describe it below (we 
give further details in Appendix [B| . 

• IV. We use the same set of fields as in approach III but focus on a different region of 
parameter space. An appropriate redefinition of the fields Hd-, Ru and Rd in approach III 
can give rise to viable down-type fermion masses while keeping the i?-symmetry intact. The 
i?-charges of L, D'^, E'^ are equal to unity, the same as in approach II and III. However, we 
redefine the three fields Hd, Ru, and Rd with R[Hd] = 0, R[Ru] = 2, R[Rd] = 2 in approach 



III in the following way: 



{Hd, Ru, Rd} -> {H'd, Hd, H'^} , (1) 



respectively, so that R[Hd] = 2, R[H'^] = 2, R[H'^] = 0. Since R[Hd] = 2 now, the coupling 
XgSHuHd is allowed and a strongly first order EWPT is possible. Furthermore, down-type 
fermion masses are generated from supersymmetry breaking similar to that in approach II 
but with an important difference. The down-type fermion masses are only suppressed by the 
masses of H'^ and H'^, not the messenger mass M^ as in approach II. Hence, viable down-type 
fermion masses without i?-symmetry breaking can be generated for high supersymmetry 
breaking scales as long as Mh^ and Mh' are only slightly above the superpartner mass 
scale. When the i?-symmetry is ultimately broken by gravitational effects leading to a 
tiny vev for Hd, there will be an additional contribution to down-type masses, which will 
generically be quite suppressed compared to the one above. 

Thus, both approach I and IV are consistent with the framework considered in this paper, and 
we will consider both possibilities. We will show, however, in section [7] that the two approaches 
can in principle lead to different predictions for EDMs. Finally, let us comment on the /z-term, 
which is allowed by the i?-charges. Here we assume that a large /i-term (compared to the EW 
scale) is absent in the microscopic theory which gives rise to the i?-symmetric theory at the 
electroweak scale. For example, string selection rules often exponentially suppress coefficients of 
terms in the superpotential which are otherwise allowed by low energy effective field theory [31] . 
Thus, it is possible to naturally have an EW scale n termj^ It is then possible to absorb an EW 
scale /i parameter by a redefinition of 5*. Although the /x-term appears in the parameters of the 
scalar potential, since we will consider the most general form of the potential consistent with 
the R symmetry, and with all mass parameters of order the electroweak scale, there is no loss of 
generality. 



^Any value of /i is technically natural. 



2.1 Theoretical Details 



We envision the following simple model of supersymmetry breaking which preserves an approxi- 
mate ^(1)_R, as already discussed in [3D]. We will, therefore, be brief and only highlight the main 
features. We assume a hidden supersymmetry breaking sector with gauge interactions including 
a U{1) gauge field W^ which gets a nonvanishing D-teim, and chiral superfields X charged under 
the hidden sector gauge group which get non-vanishing F-terms. X and W^ are assumed to have 
i?-charges 2 and 1, respectively. Supersymmetry breaking is then mediated to the visible sector 
by interactions suppressed by the mediator mass M^. Since we want to be as model-independent 
as possible, we will leave this mass scale unspecified. Although in this work we will not be 
concerned with embedding the model at the electroweak scale in a UV complete microscopic 
theory, it is not hard to envision a variety of ways, which may or may not require a very small 
value of v^jj in which this can be done. Requiring a WIMP DM candidate, however, implies 
that the supersymmetry breaking scale be large. Note that this does not automatically imply 
gravity or anomaly mediated supersymmetry breaking, as high-scale gauge mediation models 
with neutralino LSP can be constructed [32] . 

With the above considerations in mind, in addition to the kinetic terms, the Lagrangian of 
our model is given by: 



L 

Lw 

^down 



^soft 



Lw + -^down + -^soft 



j(fe (yu QU'Hu + A, SHuHd + At THuHa) + h.c. , 
jdH (yd QD'^H, + ye LE'H,) + h.c. , 
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(4) 



d^e yd 
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V2ca d^e 



m7 



W^ S, + h.c. + 



(^)^^--]-^^<^)(^^--)- 



^"^''*''(ir)«'«- 



«J 
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(6) 
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where a = U{1)y, SU{2)i, SU{3)c, ^a = {S, T, O}, and Qi runs over all matter fields. Note that 
since there are no gauge singlets with supersymmetry breaking F-terms, the standard trilinear 
parameters and the Bfi parameter vanish. Let us now comment on the operators above. 



'Still with supcrpartner masses around the electroweak scale. 



The operators in L^ include the holomorphic Yukawa couphngs for the up-quark sector. 
Operators in Ljown are responsible for down-type fermion masses. These can either arise from 
holomorphic Yukawa couplings (as in approach I) or by supersymmetry breaking (as in ap- 
proach IV described above and in Appendix [B]). The qualitative details of our analysis for the 
nature of the phase transition and WIMP DM abundance are essentially independent of the two 
approaches, hence the results obtained here apply equally well to both approaches. 

Operator (5) in Lsoft not only gives rise to Dirac gaugino masses, Mn^ = Ca m-, but also to 



mass terms for 



"he real part of the adjoint scalars Re(E(j) as well as "trilinear" couplings of the 
form (M^Sa + h.c.) ((/„ ^^ q*T°'qi) [23] ■ These new trilinear couplings will be relevant when we 
discuss the EWPT in Section |3} Operator ^ can be written for fields in real representations 
of the gauge group, hence it provides (equal and opposite) mass-squared terms for the real and 
imaginary parts of S^. Operator ([T]) generates the standard soft mass-squareds for all scalars. 
Note that the adjoint scalars S^ receive masses from all three operators above, resulting in 
different masses for their real and imaginary parts. 

Since the hypercharge adjoint is a singlet S with vanishing i?-charge, it is possible to write 
a trilinear parameter for it: 



-my 



</^'M^^l^' + h.c. 



which is, however, parametrically suppressed relative to other parameters -it is of 0{M'^^^^/Mi,) 
compared to (9(Msoft) for others. In general, there could be an anomaly-mediated contribution, 
which is also small compared to other soft terms. Higher powers of S appearing in the scalar 
potential will be suppressed even further. So we will neglect the above term (and other higher 
order powers of S) in our subsequent analysis. 

As explained earlier, a small vev for Hd is induced once a B^ term is generated. Thus, it 
is natural to have Vd ^ Vu in this framework. However, the precise value of the suppression 
depends on the mechanism of /^-breaking. Note that the usual definition of "tan /3" is not very 
appropriate, and will be discussed in Section [7} It is also important to clarify that there is no 
light (pseudo) scalar, like the i?-axion. This is because the i?-symmetry is broken explicitly by 
small amounts, and the same interactions which break the i?-symmetry (such as the B^x term) 
generate a small vev for the i?-charged Hd- Hence, the situation is similar to the MSSM with 
large tan/3 and small -B/x, in which there is no light pseudo-scalar. Finally, note that even in 
the presence of i?-symmetry breaking by Majorana gaugino masses or the B^i term, a discrete 
Z2 symmetry (call it i?-parity) is preserved making the LSP completely stable. 

As has been realized, there also exist potentially dangerous operators which must be sup- 



pressed. These include: 

'(feW'^WaY + '^.c. , (9) 

d^e —— s . (10) 

Operator ^ induces kinetic mixing between hypercharge and the hidden f/(l), and gives rise 



to an undesirable large hypercharge D-term, while operators (10) give rise to tadpole terms for 
the singlet S which could destabilize the electroweak hierarchy [33] . Suppressing the coefficients 
of operators arising from the superpotential is technically natural. Appealing to that is not 
possible for the coefficient of the operator coming from the Kahler potential, however. 

Many possible mechanisms have been discussed in order to suppress the coefficients of the 
above operators. For example, in gravity-mediated models or higher-dimensional constructions 
this may happen if S originates from a a non-abelian gauge theory (in many cases a GUT 
multiplet) which is unbroken at a scale not far from the TeV scale [301 El]- Within gauge- 
mediation models, it is possible to avoid a large tadpole term if the couplings between adjoints 
and messengers respect a GUT symmetry [26], or if there is an appropriate parity symmetry 
forbidding a large tadpole term [35] . Finally, a large tadpole can be avoided if the theory itself 
has a low cutoff A.,,, close to the TeV scale [21]. Since the focus of the paper is not on detailed 
model-building, we will be agnostic about the precise embedding of the model in a UV complete 
theory and analyze the effective theory defined by ([3])-([7]), assuming that the operators in ^- 



(10) are sufficiently suppressed, such that the singlet tadpole is of EW size. This will be crucial 



for the phase transition, as will be seen below. 

3 The U{1)r Symmetric Limit: T = Analysis 

We are now in a position to study the scalar potential of the model. With the Lagrangian given 
in Eq. ^ of Section 2.1, the scalar potential takes the form: 



= 


Vp + Vn + V^oit + ^ii , 


dW 


2 . 3 

^ a-1 



V = Vp + Vn + Knft, + V^:i , (11^ 

where 

i a=l 

Koft = mj,jH^\^ + mljH,\^ + ml\S\-' + mlT''^T'' (13) 

+ BtT^T" + tsS + BsS"^ + AsS^ + h.c , 



while Vi^J refers to the 1-loop contribution to the effective potential, which will be specified 
below. Here i runs over all fields appearing in the f/(l)ij-syninietric supersymmetric Standard 
Model, and Dy and Df are the hypercharge and SU{2)l D-terms, respectively. Compared to 
the MSSM case, the D-terms contain additional pieces associated with the SU{2)l and f/(l)y 
adjoint fields: 



Dl = ginyHu + HyHd + T^X'T) + V2 {Md, T^ + h.c.) , (14) 

Dy = g'{HiH^-HlH,) + V2{Mn,S + h.c.) , 

where r" and A" are the two and three-dimensional SU{2) generators respectively. Note that. 



as explained in Section 2.1, the D terms above give rise to new trilinear couplings in the scalar 



potential, which will be relevant when we discuss the phase transition. Also, the masses of the 



real and imaginary parts oi S = Sn + i Sj and T = Tr + iTj are split in Eq. (11). For instance, if 
Md^, Bs and Bt are real, then m|^ = ml+2B s+AM^,^, m|^ = ml-2Bs, m^^ = m^+2BT+AMl,^ 
and m^ = m\. — 2Bt- In this section and the following we will assume that there are no CP- 
violating phases in the Higgs scalar potential. We will comment on the possible presence of 
such phases, and their connection to EWBG, in Section [7| Note also that, within approach IV, 
there exist additional fields H'^ and if^, but since they are assumed to be parametrically heavier 
than the other fields (see Appendix Ib]) , their effect on the minimization of the potential can be 
neglected . 

In order to minimize the above potential, we point out two important simplifications. First, 
EW precision constraints require the triplet Higgs vev, T^ = vt, to be small, which can be 
achieved if the triplet soft breaking mass, m^, is in the multi-TeV range. Therefore, the effect 
of the triplet on the minimization of the potential must be small, and we will set wt = 0. [j In 
addition, we will analyze the potential in the /2-symmetric limit. The ?7(1)_r symmetry with 
our /^-charge assignments implies that b = Bfi = 0. Hence, if m'jj > 0, one can easily see 
that (Hd) = (i.e. we avoid spontaneous breaking of the U{l)ji symmetry), and the degrees of 
freedom in Hd effectively decouple from the minimization of the potential. Therefore, the EW 
symmetry is broken by (H^) = v, with v ~ 174 GeV, and the eaten Nambu-Goldstone bosons, 
as well as the SM-like Higgs, all arise from H^- In particular, the vacuum is automatically 
charge-preserving, exactly as in the SM. Note that when i?-symmetry is broken, there will be 
other terms in the scalar potential. However, those terms represent a small perturbation on the 
analysis based on the previous approximations, and do not affect the qualitative results. 

The zero-temperature 1-loop contribution to the effective potential, in the i?-symmetric 



^We could keep track of the small triplet vev, but setting it to zero will make the physics more transparent. 



10 



limitQis given by l^^^jp = ^ AA |i7„|^ where [371 EH] 



AA 






^^-'SUSY 
2 

mi 



y't-\{9' + 9") 



(15) 



with yt the top Yukawa couphng and rrit ^ 173 GeV the top quark mass. Note that we do not 
write contributions to AA from the A-terms, since these are forbidden by the U{1)r symmetry, 
and remain small provided the U{1)r violation is small. Thus, the size of AA is controlled by 
^susY, which can be taken as the geometric mean of the LH and RH stop masses. 
Minimizing with respect to i7°, we get the condition 

5 = mj,^+^-{g^ + g'^ + AAX)v^ + V2g'MD,Vs + Xy, = 0, (16) 

where Vg is the vev of S. We will also take advantage of the fact that Ag is expected to be 
parametrically suppressed with respect to the other parameters. So, neglecting Ag one gets the 
following relatively simple expression for Vg-. 

2tg + ^/2g'MDy 



mt 



(17) 



^5. + ^l^') 

The vacuum in the above limit thus only depends on As, nig , tg and Md^. 

The spectrum of Higgses is as follows. Due to the constraints from EW precision measure- 
ments, the heaviest of the Higgs states are expected to be the excitations of the triplet T", with 
a mass of order mrp ~ few TeV. There is also another charged Higgs and a complex neutral 
scalar, both arising from if^, with masses 

1 






<i 



^H,+ 



K + K 
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(/ + 9' 



/2\ 



;2 - V2g'MD,Vg + Xy, 



^H" 



mU + I - ^' - A^ - A| 



(19) 



m| + Xy — QAgVg, as well as a 



which could also be somewhat heavy, depending on the soft mass mjj 

Finally, there is a CP-odd (singlet) scalar with mass m^ 
pair of CP-even Higgs states that can have couplings to gauge boson pairs. The latter are linear 
superpositions of h^ = Re(if°) — v and s = Sr — Vg, and the corresponding mass matrix in the 
{/i°, s} basis reads: 

l{g^ + g'^ + 4AX)v^ + 6 



M\ 



W^g'Mo, + IXl 



[V2g'Mn, + 2X1 



mc„ + Xy + QAgVg 



4. 



(20) 



^In approach IV of Section [2J where the QD'^Hd superpotential term is forbidden, there are no terms in the 
one- loop effective potential involving Hd- In approach I, such terms are generated, but since we will always be 
in a situation where H^ is very small, these terms play only a minor role in the minimization of the potential. 
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where we included the dependence on As for completeness, although we expect it to be negligible. 



We also included in the (1,1) entry the contribution from 5 in Eq. (16): although it vanishes 



in the zero-temperature vacuum, it will give a non- vanishing contribution at finite temperature, 
when V and Vg are not identified with their zero-temperature values, and should be included 



in the analysis of the next section. Note that the CP-even states in (20) depend on the same 



set of microscopic parameters that determine the vacuum, while the remaining Higgses depend 
on additional parameters that do not affect the vacuum (in the limit defined above). Thus, the 



stability along the excitations of Eqs. (18) and (19), as well as the CP-odd singlet and the triplet 
directions, is easily guaranteed by the choice of mjj- , m| and m^, respectively. 



The neutralinos and charginos also play an important role in the collider phenomenology, 
the DM question, and the EWPT. Since the neutralinos are Dirac in nature in the i?-symmetric 
limit, the mass matrix in the basis {iB, iW^, H^ ; T, S, H^} takes the form: 



M^o 







X 



N 



Xn 




where 



X 



N 




(21) 



and we again set the triplet vev to zero. 

Similarly, for the charginos we have in the {T~ ,iW~ ,H^ ; T^,iW^,H^} basis: 





M^± 



X, 



T 



Xc 




where 



Xr 




(22) 



At zero-temperature Vs is given by Eq. (17). Note that the neutralino and chargino mass matrices 



depend on two parameters that do not enter in the determination of the vacuum: Ay and M^^. 
Furthermore, it is easy to find that the three chargino mass eigenvalues are given by 



^xt 



M, 



D2 5 



X2,3 



2 , Ml,, + (/ + 2\iy + Xlvl 



± J [Ml^ + {g^ + 2\l)v^ + Xlvl] ' - 4 (a.M^.t;, + ^/2g\TV^^ ' 



(23) 



(24) 
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We see from the first equation tliat we must require M^ij > 104 GeV to avoid a chargino 
lighter than the direct LEP chargino bound [M]. However, we also note that we cannot take 
M^2 arbitrarily large, since in that limit one finds m ± ^ XsVs[^ + C>{v/M£)-^,Vs/M£,-^)], which 
is typically smaller than the chargino LEP bound. Thus, the chargino bound selects a region 
in Md2 which, together with the region in Md^ and A^ preferred by the minimization of the 
potential (and the requirement of a strongly first-order phase transition to be discussed in the 



next section) leads to a particular composition for the lightest neutralino, from Eq. (21). In 
particular, one finds a sizable Higgsino component, which has interesting consequences for the 
computation of the DM relic density and for direct DM searches, to be described in Section [5] 

4 A Strongly First-Order EWPT 

We now turn to the finite-temperature potential, focusing in particular on highlighting the 
reasons that lead to a strongly first-order electroweak phase transition (EWPT). We assume for 
concreteness that the early Universe had a high-enough temperature. The temperature decreases 
with the expansion of the Universe, and EWPT takes place at temperatures of 0(Tew) = Mew, 
which is the object of our study in the following. The main effects arise already at "tree-level", 
so that it will be sufficient to simply consider the thermal masses induced by the plasma. A 
more refined analysis based on the full finite-temperature effective potential is possible, but is 
not expected to change the conclusions other than in small details. This is to be contrasted 
with the case of the MSSM (or the SM), where carefully considering up to two-loop ("daisy" 
improved) effects [lO] is essential to determine the strength of the phase transition (see [H] for 
a recent review). 

We continue -as in the previous section- to consider the limit where the triplet Higgses are 
heavy so that the triplet vev is sufficiently suppressed to be consistent with the constraints on 
the Peskin-Takeuchi T-parameter. Therefore, these particles are effectively decoupled from the 
plasma. However, we will assume that the four real degrees of freedom in H^ and the CP-odd 
singlet scalar Sj can be considered as light degrees of freedom at the temperatures of interest (all 
the masses below ~ ttTc, where Tc is the critical temperature). The two CP-even Higgs states 
from the h^-s system, and the neutralinos/charginos are also light, hence active in the plasma 
during the EWPT. 

In the f/(l)ij symmetric limit considered in the previous section, the minimization of the 
potential at finite temperature effectively involves only two degrees of freedom: (p = H^ and 
the singlet scalar, 0^ = Sr (at zero temperature, (0) = v and (0^) = Vg). The thermal masses 
can be read directly from the terms involving (i/°) = v and (Sr) = Vg given in the previous 
section. We use the high-temperature expansion for the effective potential, so that the leading 
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T-dependent terms are obtained from 

1 



Vt 



24 



T' 



2m\jo + 2m^± + m^ + TiMjj + 6m^ + 3m| 



48 



2Tr(>l],oA^^o) + 2Tr(A^LA^^±) + 12m 



■^2 j2 



= C(i,T (f) + C52T 0^ + C5T 0^ + const. , 



where we can drop the field-independent constants, and 



32 y 32 y 4 ^* 8 4 ' 3 ^ ' 



A? , c.s = ^ ^'Mb, 



4^2 



(25) 

(26) 
(27) 



We also included the contributions from the electroweak gauge bosons and the top quark: m^ = 
|f7^0^, ml = ^{g'^ + fi''^)0^ and m^ = yt(f) (in the tan/3 = cxd limit). The finite-temperature 
potential then takes the form 

V = {m^ct)^ + A0^) + (2t,0, + m^02 ^ ^^^3) ^ (2S,0,02 ^ ^^^2^2^ ^ (28) 

where the effective coefficients (denoted by a tilde) are given in terms of the parameters of 



Eqs. (11) and (26)-(27) by 



~ 2 
m 



m: 






as 



ts + \csT^ 



A 



A. 



1 

8 

A? 



{g^ + g'^) + i AA 



(29) 



The potential of Eq. ( 28 ) captures the physics to a good approximation even when allowing for 



a non-zero but small triplet vev, and for small U{1)r violating effects that can induce a small 
Hd vev. 



We have split the terms in Eq. (28) into terms that depend only on 0, terms that depend 
only on (ps, and terms that mix and (ps- In order to understand the minimization of this 
system, and to highlight the underlying mechanism that leads to a strongly first-order phase 
transition, it is useful to start from the situation where the mixing terms vanish. In addition, 
we will make use of the fact that A^ is expected to be small [see comments following Eq. (Is])], 
and we will neglect it in the following. Then one can see that the potential defines two scales: 
— m^/(2A) (= 0^ if m^ < 0), and 0s = — ts/m^. We are interested in understanding the effects 
of a non- vanishing singlet vev on once the mixing terms are taken into account. This can be 
understood in a transparent way by writing an effective potential for after "integrating out" 
0s at tree-level (i.e. using the 0^ EOM). The mixing terms introduce a further scale into the 
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problem, 0^j, = rhl/Xg, that defines a transition field value for 0, below which the mixing terms 
are a small contribution to the potential compared to the pure singlet terms, and above which 
they dominate over the pure singlet terms. For instance, in the region where ^ 0tr one has 

1 + I ^^^ - 1 1 ^ + o' ""' ' 



t. 



ml 



Aoto 



^ 



and replacing back in Eq. (28) one finds the effective potential for 

?2 



v;ff(0) 



-A. 

ml 



+ m^g02 + X 



offS 



+ w 



6a 



where 



m 



eff 



m^ 



2osts A^t^ 



?Ti? 



m: 



(30) 



for 6 ^ 



^tr , 



(31) 



(32) 



A 



eff 



A 



m? 



Xsis - asml 



(33) 



We see that one important effect of the singlet vev is that it gives a negative contribution to 
the effective quartic coupling near the origin. In fact, as we will see, the condition Aeff < 
characterizes the potentials that exhibit a (strongly) first-order phase transition at tree-leveln 
The "instability" signaled by the negative effective quartic coupling indicates that there is a 



minimum far away from the origin. In the small field expansion of Eq. (31), this "far away' 



minimum can be understood as arising from balancing the quartic operator against the tower of 



higher-dimension operators in Eq. (31). 



That the "far away" minimum exists in the above situation is guaranteed by the fact that 
(softly broken) supersymmetric potentials are bounded from below. One can see this explicitly 
by considering the case where ^ 0tr- In this region, the mixing terms in Eq. (28 ) dominate over 



those that depend only on the singlet, and therefore the singlet vev is now given approximately 

by 



A. 



l + O 



so that 



v;ff(0) 



const. 



?Ti^ 



as 



+ A0^ + O(l/02 



for (h ^ 



Jtr 



(34) 



(35) 



Thus, the original (positive) quartic operator controls the large-field behavior of the potential, 
while the singlet vev affects only the quadratic and other subdominant terms. 



'This point has also been noted in 
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The above arguments suggest that there is a first-order phase transition, and that it can 
easily be strong. Indeed, if m^g > (either at T = or due to the thermal mass contribution), 
then there is a local minimum at the origin, that can be separated by a barrier from the "far 
away" minimum described above. For this to be the case, one must also have that the effective 
mass, ■nigfj, should not be excessively positive, or else it will overwhelm the effect of the negative 
quartic operator. On the other hand, if this squared mass is negative at T = 0, the thermal 
mass contribution will make it eventually positive, and there will exist a barrier provided only 
that Aeff < 0. If the symmetry breaking minimum is the global minimum, then the main effect of 
increasing the temperature is to raise the global minimum to be degenerate with the minimum 
at the origin at some critical temperature Tc. It is worth emphasizing that this mechanism does 
not rely on the existence of a term cubic in 0, and is therefore qualitatively different from the 
way the phase transition works in the MSSM (or, had the Higgs been light enough, in the SM). 
As we will see, this allows more easily for the phase transition to be strongly first-order. 



Potentials of the form of Eq. (28) were analyzed in Ref. [19], in the limit that the temperature 
dependence of rh^ and ig is neglected, i.e. setting m^ = m| and ig = tg in Eq. (29). The 



model studied in that work was different from ours, and therefore the relation of the effective 
parameters to the underlying ones is different in the two cases. However, at the level of the 
effective parameters the two scenarios are identical, and we simply summarize the main results 
of the analysis in [19]. Defining the critical temperature Tc and the critical vev (pc (= (0) at Tc) 
by: 

(36) 
(37) 



one finds 



where c^f, was defined in Eq. (26), and 



y(0, 


„T,) = y(0 = O,T,), 


dV 

d(j) 


= 0, 


1 / , 1 


/\gLg 

msj^ag 


C0 


26), and 


F(02) 


= -2as(j)s - 


A,02_2A02 . 



(38) 
(39) 



Here, it is understood that 



2ts + ^/2g'MD^(P 
2(m|^ + A202) 
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which should be compared to Eq. (17). Based on the requirement that 0^ > 0, the authors of 
Ref. 



derived a necessary condition for a first-order phase transition that can be seen to be 
equivalent to the requirement that Xcs(T = 0) < 0, where Agfr is the effective quartic coupling in 



the small expansion, as defined in Eq. (33). Our previous comments then give a transparent 



interpretation for this condition in regards to the existence of a barrier separating a minimum 
at the origin from a "far away" minimum. Furthermore, it clarifies why this is only a necessary 
-but not a sufficient- condition for obtaining two degenerate minima separated by a barrier, and 
allows a generalization to the case that the temperature dependence in m^ and tg is included. 
As already mentioned, if the zero-temperature effective mass parameter, mlg is positive and 
excessively large, it may overwhelm the negative quartic operator in the small region. What 
is required is that there exist an "intermediate" region where the quartic operator can become 
relevant: at very small 0, the potential is determined by the quadratic term, while at large the 



higher-dimension operators [in the language of Eq. (31)] dominate. For instance, the coefficient 
of the 0^ term in Eq. (31) is given by Geflf = As(asm^ — Xsts)^/rh^ > 0. Thus, we see that the 
"far away" minimum is of order 0^ ~ a/— Aeg/Geff, which is obtained by balancing the 0^ and 
0^ operators (the actual value of the "far away" vev depends on the full "tower of operators" in 



Eq. (31), but the above approximation captures the qualitative dependence up to an order one 



coefficient). Therefore, by requiring that the quadratic term not be larger than the quartic one 
for 0c ~ a/— Aeff/Geff, one finds that one must also have GcsttiIs ~ -^cfi- Strictly speaking, the 
above relations (Aes < and a not too large m-gg) should hold at the critical temperature, where 
the existence of the barrier is crucial. However, an approximate criterion for a first order phase 
transition in terms of microscopic parameters is 



A 



efflr=0 



A 



TTl 



Sr 



A.ot.o 



a^nii 



G^effW^efflr^o 



As(as"^|^ - Xsts 



m' 



Sr 



rn^ — 



<0 



ZQjcZq 



+ 



m 



Sr 



Kt'g 



m 



Sr 



Note that if w^efflT=o ^ ^^ ^^^ second condition is automatic. 



< 



X2 I 
^eff I T=0 



(40) 



The condition that T^ > [as obtained from Eq. (39)] was also considered in Ref. |19j . 
This condition also has a simple interpretation: it is equivalent to requiring that the "far away" 
minimum be the global one (at T = 0). The EW-preserving minimum has = and 0^ = 



—ts/m"^ , which leads to an associated potential energy Vq = —t1/m?g [see Eq. (31)]. At the 
"far away" minimum one has 



V{v,Vs 



Vrs 



[l + XgV^ jvfi^ 



Acff + 2AA, 1 + 



Sr} 






m 



Sr 



(41) 
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Figure 1: The hatched region in the tg-X^ plane leads to a strongly first-order phase transition (red) and is 
allowed by the Higgs LEP bound (green). We have fixed Md, = 35 GeV, msj^ = 100 GeV and Msusy = 2 TeV. 



where Aeg, as given in Eq. (33), is evaluated at T = and Vs is given in Eq. (17). Requiring 



that this expression be negative (to ensure that the EWSB minimum is the global one), results 
in the condition 



A < 



(As^s 



a,m|^)2 



ml. 



[ml. 



+ Kv^Y 



(42) 



We show in Fig. [T] how the various constraints determine the physically interesting region of 
parameter space. As emphasized before, the vacuum structure depends only on four microscopic 
parameters: As, ""^1^, tg and Mr,^ (if As can be neglected). We show a projection on the tg- 
Xs plane (note that we use the microscopic As, not the effective As), for fixed M^j = 35 GeV 
and m| = (100 GeV)^. In this figure we neglect the temperature dependence in ml and ts, 



just as was done in Ref. [19], but we will make no approximations in the numerical studies of 
subsequent sections. The light blue region (to the right of the boundary marked as "0c exists") 
corresponds to the condition 0^ > 0, which coincides with the condition AegL^Q < 0. The 



dashed line corresponds to the boundary where Eq. (42) is saturated, so that to the right of this 



line the theory does not break the EW symmetry. The interesting region then lies between the 
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Figure 2: Left panel: lines of constant m^jj/u^, as defined by Eq. (32). Right panel: lines of constant AoS, as 



defined by Eq. (33 1. Both quantities are evaluated at T = 0. The shaded (light blue) region corresponds to that 



consistent with a strongly first-order EWPT and a heavy enough Higgs mass, for the same choice of parameters 
of Fig. [1] 



two previous boundaries. In addition, we denote in a red shade the region where the EWPT is 
strongly first-order, as defined by the criterion (pc/Tc > l/v^ (recall that our normalization is 
such that at T = 0, (0) = 174 GeV). It is clear from the figure that the condition Xes\T=o ^ ^ 
leads almost automatically to a strongly first-order EWPT. This point has been emphasized 
more generically in |12]. We have also indicated in green the region where the lightest CP-even 
Higgs has a mass rrih > 114 GeV, where 



m. 



TrA^2 



H 



{TrM%y - 4:DetMjj 



(43) 



and Aiff is given in Eq. (20) (with 6 = 0). Although this condition is overly restrictive, since the 



lightest CP-even Higgs state has in general a singlet component, it illustrates that the Higgs LEP 
bound is easily consistent with a strongly first-order EWPT. It should be noted, however, that 
obtaining a heavy enough Higgs mass requires that the 1-loop radiative corrections be significant. 
In particular, in the tan/3 = oo limit, the presence of the singlet does not affect the Higgs mass 



compared to the MSSM [see the expression for A in Eq. (29)]. However, it is also important that 
there is no need for either top scalar superpartner to be light in order to achieve a strongly first- 
order EWPT, unlike in the light stop scenario of the MSSM. Hence, both the LH and RH stops 



can be taken somewhat heavy. In the figure, we have chosen mi 



mz 



MsusY = 2 TeV. The 



hatched region corresponds to models that exhibit a strongly first-order EWPT and are allowed 
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by the LEP Higgs bounds. The shape of the allowed region is similar for other choices of M£,^ and 
m|^. However, as we will see, there are further constraints coming from the neutralino/chargino 
sector (in connection to DM) that prefer relatively small values for M^,^. Note also that A^ is 
required to be sizable although, as we will see, DM constraints have a preference for the lower 
range of values. 

In Fig. [2]we show again the allowed region, for the same parameters of Fig. [l} and superimpose 
the lines of constant ^«eff/'^^lT=o (^^^^ panel) and \eE\T=o (right panel). This shows that the 
allowed region is indeed characterized by Acfrl^i^Q < 0, as well as moderate ^efflT=o ^^ ^^^ sense 
that the effective mass parameter is at most of order v"^ (it is also clear that this does not 
represent an overly restrictive condition). It can also be checked that the allowed region satisfies 
Geff'Tigfj/A^gl^.^Q < 0.9, thus validating the physical picture described in the paragraph before 



Eq (40). 



5 Pseudo-Dirac Dark Matter - Constraints and Signals 

We now discuss various aspects of Dark Matter physics in detail. We assume that the scale of 
supersymmetry breaking is high enough that the lightest i?-parity odd superpartner, which is 



a DM candidate, is a WIMP -usually the lightest neutralino, ^° This is in contrast to models 
with low-scale supersymmetry breaking where the gravitino is likely the DM candidate. Thus, 
the relic abundance of the DM candidate is determined by a freeze-out calculation]^ and the 
LSP is also expected to be seen in direct-detection experiments. In the U{1)r symmetric limit, 
the neutralino is Dirac in nature and scatters with nuclei either via s-channel squark exchange 
if it is gaugino like, or via t-channel Z-exchange if it has an appreciable Higgsino component, 
through a vector-vector interaction which is not suppressed by the DM velocity. As such, these 
contributions are subject to extremely stringent constraints. Hence, a Dirac neutralino DM 
candidate is essentially ruled out. 

However, in the presence of (small) /?- violating operators which are necessarily present, 
the two degenerate mass eigenstates of the Dirac DM split into two Majorana fermions with 
(slightly) non-degenerate masses. Hence the DM in this framework is pseudo-Dirac in na- 
ture. This can naturally allow for consistency with the latest direct-detection constraints from 
XENONIOO |13] if the splitting between the eigenstates is larger than the typical momentum 
transfer (of O(IOO) keV) in a direct-detection experiment, such that the transition from one 
state to the other is kinematically forbidden. In our framework, we assume that the /^-breaking 

^''Note that i?-parity is still conserved in the model even when C/(1)_r is broken by Majorana gaugino masses 
and/or the B^ term. 

^^Assuming a standard radiation-dominated era after the end of inflation until after BBN. 
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is such that the two Majorana eigenstates are spht by 5m > 100 keV (typically a few GeV, which 
could arise from anomaly mediation for example). Thus, the pseudo-Dirac DM is Majorana-like 
for the purposes of direct-detection. Note that this is quite different from the scenario of inelas- 
tic DM (iDM) IH], where the two states are split by an amount comparable to the momentum 
transfer in direct detection experiments. The situation regarding relic abundance can be com- 
pletely different, however, since the freeze-out temperature Tp of the LSP (of (9(1 — 10) GeV 
for a WIMP) is much larger than the typical momentum transfer in direct-detection. Thus, as 
long as 6m < Tp, which can naturally occur in our framework, the DM properties are similar to 
those of a Dirac particle in regards to the relic abundance, as a result of coannihilation processes 
being in equilibrium. 

5.1 Relic Abundance 

As mentioned earlier, we assume for concreteness that the universe passes through a "standard" 
radiation dominated era beginning from a high temperature due to reheating after inflation, 
which persists until the time of matter-radiation equality. The relic abundance of a WIMP DM 
candidate is then determined by a thermal freeze-out calculation. The annihilation channels 
relevant for computing the relic abundance of DM after freeze-out can be divided broadly into 
three classes: a) annihilation into fermions, b) annihilation into Higgs/W/Z states, and c) 
annihilation into gluon/photon final states. 

What can be said about the above channels within the class of models considered? The 
region of parameter space consistent with a strongly first-order EWPT and LEP bounds for 
Higgs and chargino masses (as well as the direct detection constraints to be discussed in the 
next subsection), typically gives rise to a DM candidate with a mass such that the annihilation 
channels in (6) above are kinematically forbidden or suppressed. Furthermore, since DM is 
supposed to be neutral with respect to electric charge and color, DM annihilation channels in 
(c) occur via loop effects which may become important only when the annihilation channels in 
(a) are suppressed, such as due to helicity suppression for s-wave annihilation of Majorana DM. 

Therefore, for the pseudo-Dirac DM with 6m < Tp, arising in the scenarios considered in 
this work, only s-wave annihilation into fermion final states is relevant, which will be assumed 
from now on. This class of annihilation channels can be further subdivided into three types 
depending upon the particle exchanged in the s or t channel. Thus, the annihilation to fermion 
final states (//) may proceed through squark/slepton exchange, Z exchange, or Higgs exchange. 
Since the squark/slepton exchange contribution depends on their (model-dependent) masses and 
does not qualitatively affect any other aspect of the physics considered in this paper, we assume 
that the masses of squarks/sleptons are heavy enough such that their contribution to the relic 
abundance is subdominant and can be neglected. Furthermore, the Higgs exchange contribution 
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is suppressed by the masses of the ferinionsj^ hence the dominant s-wave annihilation mode is 
through Z-exchange. The Z-exchange contribution is really a co-annihilation contribution since 
Xi and X2 (the two semi-degenerate Majorana states forming the pseudo-Dirac DM candidate) 
co-annihilate through a Z into fermion pairs in an s-wave process. Also, since only the Higgsino 
components of Xi and X2 couple to the Z, the relic abundance of the LSP is correlated with 
its Higgsino component. This can be seen more concretely by considering the pure Dirac {R- 
symmetric) limit. The coupling of Z to ipi (the lightest Dirac neutralino in four- component 
notation) is given by: 



V9^ + 9'^ 1 2^ ^.YiciPL + crPr) iJi , (44) 

where cl = \U^ fr P and cr = \U^ f^ p, with Ul and Ur the unitary matrices diagonalizing the 



Dirac neutralino mass matrix, Eq. (21), i.e [/lXn Ur = diag(m^.). Thus, cl and cr correspond 



to the Hd and Hu components of ipi, respectively. For a pure Dirac neutralino, the thermally 



averaged vector- vector annihilation cross-section to fermions arises from the vector piece of (44), 
and is given in the non-relativistic limit by: 



where we defined 






(46) 



The sum runs over all SM (Weyl) fermions lighter than ipi (typically including the bottom, but 
not the top quark), and gj. = T|. — 2s'^Qf^ is the standard coupling of fi to the Z gauge boson. 
We neglected phase space suppression factors that are very small even for the bottom quark, for 
the typical neutralino masses of interest. This expression is modified slightly when i?-breaking 
is included, so that the process now corresponds to co-annihilation of x\ and x\i the masses of 
which are split by 5m. In the limit that 5m <^ m,^^, we can write approximately 

{^v)f}:^Jf ^ {av)r^;:Jf^{2 + 5f{l + 5f/'e-^-', (47) 

yoff 

2 + 2{l + 5fl^ e-^- ' ■ 5 ^"^ 



9cs 



m,^o 



Here xp = m,^o/Tp with m^o ^ m^^ and Tp the freeze-out temperature of Xi and X2) which has 
to be determined self-consistently. 
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Figure 3: Contours of the LSP relic-abundance rj^o/i^ in the tg-Xs plane for the following choice of parameters. 
Left panel: Md, = 35 GeV, M^, = -110 GeV, ms^'= 100 GeV, MgusY = 2 TeV; Mi = 5 GeV, M2 = 10 GeV, 
b = Bfi^ (40GeV)2. Right Panel: Md, = 60 GeV, Md^ = -110 GeV, tosh = 100 GeV, MgusY = 2 TeV; Afi 
= 10 GeV, M2 = 20 GeV, B/j, = (40GeV)2. The shaded (light blue) region in the plots corresponds to that 
consistent with a strongly first-order EWPT and a heavy enough Higgs mass (in the i?-synimetric limit, which 
should not change much if the small i?- violation due to b is also included in the minimization of the potential) . 



Fig. |3]shows contours of the LSP relic-abundance for two choices of parameters (see also Fig. |5] 
for further motivation for these choices). The LSP relic abundance in the left panel is always 
a small fraction (< 1%) of the total DM abundance determined by WMAP [15]. However, the 
LSP relic abundance in the right panel is typically an 0{1) fraction of the total DM abundance, 
and could even account for the entire relic abundance in a region of parameter space. The 
main reasons for the different relic abundances in the two plots above is the following. First, 
the LSP masses for the choice of parameters in the left panel are close to m^/2, as seen from 
the left panel of Fig. [4|P^ causing a resonant enhancement in the annihilation cross-section. No 
such resonance enhancement is present in the right panel since the LSP masses are larger and 
essentially outside the "resonance" region (see right panel of Fig. |4]). Also, the quantity S in 
the left-panel is generically smaller than that in the right panel. Although, the precise value 
oi 6 = ("^x° ~ rn^o)/ 171^0 varies throughout the plots, 6 is roughly given by Mi/m^o, which is 
~ 1/9 in the left panel and ~ 1/6 in the right panel above. Since the relic abundance depends 
exponentially on 5 through Eq. (|47|), and xp ~ 25, this also has a non-trivial effect on the relic 



^^The DM is typically not heavy enough to annihilate into tt. 

^•^Fig. H also illustrates that the mass of the lightest neutralino (LSP) for these two choices of parameters is 
consistent with the LEP bound on the invisible Z-width. The lightest chargino mass (> 105 GeV) for these 
choices of parameters is also above the direct LEP bound. 
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Figure 4: Contours of the lightest neutrahno (DM candidate) mass in the tg-Xg plane for the two choices of 
parameters of Fig. ^ The shaded (light blue) region corresponds to that consistent with a strongly first-order 
EWPT and a heavy enough Higgs mass, in the i?-symmetric limit. 



abundance. Hence, demanding that the LSP abundance accounts for all the DM in the Universe 
favors a relatively "large" value of 6 , implying a relatively "large" amount of /^-breaking. 



However, since the i?-breaking is technically natural^ all effects arising from the breaking of 
the i?-symmetry are naturally suppressed by powers of 6 and loop factors. Furthermore, R- 
breaking by Majorana gaugino masses does not affect the scalar potential and hence the shape 
of the shaded (light blue) region in the figures at leading order. 



5.2 Direct Detection 

Having established that the LSP relic density can be a sizable fraction of the energy content of 
the universe, we now turn to its detection. Since the pseudo-Dirac LSP behaves like a Majorana 
particle for direct detection, the two main channels through which it can scatter off a nucleon 
is through an axial-vector interaction or a scalar-scalar interaction. The axial-vector interaction 
goes like v"^, where v is the present velocity of DM in the surrounding halo (~ 10^^ c); hence, 
the dominant interaction is via a scalar-scalar interaction, through Higgs exchange. The Higgs 
exchange contribution arises only if the lightest neutralino has a non-trivial Higgsino component. 
What is the typical size of the Higgsino component in our framework? The parameter space 
for a non-trivial Higgsino component is correlated with that giving rise to a strongly first-order 



^■^In the sense that one has an enhanced symmetry in the limit when the coefficients of the i?-breaking operators 
vanish. 
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Figure 5: Contours of the spin-independent direct detection cross-section (in units of 10 ''^cm^) in the tg-Xs 
plane for the following choice of parameters. Left panel: Mdi — 35 GeV, Mjj^ = -110 GeV, mg^ — 100 GeV, 
MsusY = 2 TeV; Ah = 5 GeV, Ah = 10 GeV, Bfi ^ (40GeV)2. Right Panel: Md, = 60 GeV, A/d, = -110 
GeV, ms^ = 100 GeV, MgusY = 2 TeV; Mi = 10 GeV, Ah = 20 GeV, Bfi = (40GeV)2. The shaded (light 
blue) region in the plots corresponds to that consistent with a strongly first-order EWPT and a heavy enough 
Higgs mass, in the i?-symmetric limit. 



EWPT in the model. This can be understood as follows. Fig. [T] shows that the existence of a 
strongly first-order EWPT with a heavy enough Higgs mass to satisfy the LEP bounds requires 
As ^ 0.6. However, it turns out that the coupling of Xi to h^, which is the Higgs state with SM 
couplings to nucleons, increases as one increases A^. Since the Xi-^u coupling is given by 



Uy.0 hO 



9 

V2 



^w + ^^b) U^. + i^sUs + XTUf)U^^, 



(4^ 



where Ua corresponds to the "a" component of Xi, ^x^\hP is linearly related to the Higgsino 
components of the LSP, i.e. Ufj and Uj^ . Therefore, the lower bound on As for a strongly 
first-order EWPT places a lower bound on the Higgsino components of DM, thereby implying a 
lower bound on the Higgs-exchange spin-independent direct-detection cross-section, cxsi. 

Taking into account the exchange of the two CP-even Higgses with couplings to fermions, 
the spin-independent nucleon/DM elastic scattering cross-section can be written, in the non- 
relativist ic limit, as 






^^X?,/i°^hNN 



TT 






\^H,hu 



mi 



^ih ^ 



2 \2 
H) 



m%ml 



(49) 



where h and H are the CP-even Higgs mass eigenstates that have an h^ component [see Eq. (20)] 
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Here U denotes the unitary matrix that diagonahzes the {h^, s} system and U^ f^^, U^^^^^ denote 
the h^ content in h and H, respectively. We also have that the effective Higgs-nucleon coupling 
is ghNN ~ 0.3((7m7v)/(2mi4/), with m^ the nucleon mass [l6]. In Fig. [s] , we show the contours 
of constant asi = (J% for the two different choices of parameters for which we showed the 

Xi 

relic density in the previous section. Since the Majorana nature of DM is crucial for direct 
detection bounds, the set of parameters includes i?-breaking Majorana gaugino masses and the 
B/i parameter (these were also included in the relic density contours). The shaded (light blue) 
region is, however, drawn in the i?-symmetric limit for simplicity, as the small i?-breaking has a 
minimal effect on the size and shape of this region. 

Fig. [5] shows that a subset of the shaded (light blue) region is consistent with the latest 
XENONIOO upper bound (7 x 10~^^ cm^ for Mdm ~ 50 GeV) for both choices of parameters. 
The allowed region is smaller in the right panel of the figure compared to the left one, mainly 
because Uu is larger for the choice of parameters in the right panel compared to that for the 
left, and because the lightest Higgs mass is approximately equal in both plots. However, it is 
worth noting that due to the lower bound on the direct-detection cross-section in the model, as 
illustrated in the plots, it is expected that the next round of results from XENONIOO should 
be sensitive to this class of models, especially if the LSP accounts for an 0{1) fraction of DM. 
If, on the other hand, the LSP accounts for a negligible fraction of DM, then the detectability 
of the signal depends on the product p{;^cai "^si since it is this combination which determines the 
rate at a direct-detection experiments^ 

Discussion: It is important to properly interpret the plots in Figs. |3] and [5] showing the results 
for relic-abundance and direct-detection. The point to note and appreciate is that depending 
on where one sits in the parameter space, the relic abundance and the expected signal for direct 
detection can vary considerably. One possibility, as seen from the plots in the right panels of 
Figs. [3] and |5| is that the LSP of the model accounts for the entire DM relic abundance and 
also provides a signal in the next round of direct-detection experiments. However, it is perfectly 
possible that the LSP abundance is a small fraction of the DM relic abundance as seen in the left 
panel of Figs. |3j This would imply that the dominant component of DM has a different origin. 
However, a direct-detection signal is nevertheless possible in this case if asi is large enough 
(greater than the XENONIOO bound of 7 x 10"^^ cm^ for moM = 50 GeV) and the local LSP 
abundance in the halo is appreciable, since the experiments are only sensitive to the product 
pfomi^si- Fig. [5] shows that large SI cross-sections are possible. Similarly, an indirect-detection 
signal in the form of cosmic-ray neutrinos from the sun for LSPs with a small relic abundance 
may be possible (see next subsection). 

The LSPs could also provide an 0{1) fraction of DM, so that DM consists of more than one 



15„LSP 



Piocai i^ ^^^ density of the LSP at the Earth's position. 
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non-negligible components. For example, axions provide an elegant solution to the strong CP- 
problem and also provide a good candidate for DM. They also naturally arise within string theory. 
Hence, it is a natural possibility that both WIMPs and axions provide an 0{1) contribution to 
the total DM in the Universe. In fact, it has been shown that this can occur in well-motivated 
particle physics frameworks arising from string theory [171 H8] . Finally, the plot in the right panel 
of Fig. [3] shows that it is also possible for the LSP abundance to be larger than the total relic 
abundance. However, the implicit assumption in these plots is that the contribution from the 
annihilation channels mediated by t-channel squark/slepton exchange is negligible. The above 
assumption is relaxed if the squarks/sleptons are light enough, so that these regions may also 
become viable in that situation (as we will discuss in Section [tI light sfermions can be consistent 
with EWBG and bounds on EDMs). 

Thus, even if the above framework provides a correct description of electroweak scale physics, 
the precise region of parameter space selected by Nature can only be determined by detailed 
measurements from a combination of experiments and observations in particle physics and as- 
trophysics. 

5.3 Indirect Detection 

In addition to "directly" detecting the WIMP DM candidate by its recoil when it scatters with 
nuclei in a detector, another possible observable signal of WIMPs arises from the production 
of electrons, positrons, antiprotons, photons and neutrinos from LSP annihilation inside the 
Galactic Halo. This "indirect" detection of DM has received a lot of attention in recent years 
following the resuhs of PAMELA ^ and FERMI [SD]. 

What are the prospects for indirect detection within our framework? The WIMP DM candi- 
date in the allowed region of parameter space tends to be light, i.e m^o < rriw Also, although 
during freeze-out, Xi ^-cts like a (pseudo) Dirac particlej^ it acts like a Majorana particle at 
present since the kinetic energy of Xi is not enough to overcome the mass difference between 
Xi and X2- Therefore, the dominant co-annihilation contribution to the cross-section giving rise 
to the relic-abundance in Fig. [3] is not operative. It follows that the signal for cosmic ray elec- 
trons and positrons is significantly below that seen by PAMELA or FERMI. The PAMELA and 
FERMI signal must have a purely astrophysical interpretation within this framework. A similar 
statement holds for cosmic 7-rays produced by Xi Xi annihilations. The DM annihilation signal 
is too weak to explain the 7-rays observed by FERMI [5T|, hence those should have a purely 
astrophysical explanation as well. 



^^In the sense that kinetic energy of Xi during the freeze-out era is large enough to allow xi and X2 ^o 
CO- annihilate. 
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What about cosmic ray neutrinos? Here, the situation is quahtatively different for a number 
of reasons. As the solar system moves in the galactic halo, the WIMPs occasionally scatter with 
the nuclei in the sun and can lose momentum and become gravitationally bound with it. For 
a wide-range of choices of interaction parameters, an equilibrium can be established between 
the annihilation and capture rate over the lifetime of the sun. Among the various annihilation 
products of WIMPs inside the sun, neutrinos are unique since only they can travel to the earth 
without significant absorption. Another important difference between indirect detection via 
neutrinos and via other cosmic rays such as photons or electrons is that the signal depends not 
on the annihilation cross-section {crv), but on the spin-independent (SI) and spin-dependent 
(SD) scattering cross-sections with nuclei (H and He in the sun)j^ 

Since the neutrino signal depends on both SI and SD cross-sections and there are very 
stringent bounds on the SI cross-section from XENONIOO, having a large SD cross-section is 
crucial for a detectable signal. The effective operator (in four-component notation) and SD 
scattering cross-section against nucleons arising from t-channel exchange of Z bosons is [52] : 

4? = K(xh'l'xl){qi,7'q), (50) 

K = -^S^nf ."^^ , ) , Vn = \U^f-\U^J\ 

\^^W'^W"''Z / 
94 171% m^o 

-^1 IT [mN + m^o)^ 

where using the results in p^, one finds approximately a^ ^ 0.705 Vu, an ~ —0.555 Vii, leading 
to a^^ ~ 1.77 X l{i~^'^ Viicm^ . Since x\ iii the allowed region has a non-negligible Higgsino 
component which couples to the Z, Vu is an 0{1) number giving rise to a^Jp > 10~^^cm^. 
The typical SD cross-section in the allowed region is much larger than the SI cross-section 
{a% < 10""''^ cm^) because the latter is dominated by Higgs exchange. Relative to the SD-cross 
section, the SI contribution is suppressed by two powers of the small effective coupling of the 
Higgs to the nucleon, ^fhNN ~ lO^'^. Since the LSP mass within the framework is expected to be 
< mwi the relevant constraints come from Super-K [53], which can be satisfied. Also, within the 
framework the channels leading to copious neutrino production (VT+W^^, t^t~ , tt) [5l] are not 
expected to be significant for the same reason. Following [55J, and the projected sensitivity of 
the IceCube/DeepCore experiment [5B], however, it is expected that this class of models should 
be observable in the near future in neutrino telescopes, if the LSP forms an (9(1) fraction of DM. 
As for direct detection, if the LSP is a negligible fraction of DM, then its detectability crucially 
depends on the combination which appears in the capture rate: {Pio<^[Cr^? rT)/m'^o- 



17i 



This is true if the capture and annihilation rates have reached equihbrium, which holds true in our models. 
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Until now we have focused primarily on the LEP and relic-abundance constraints, as well as 
on the signals in direct and indirect detection experiments. We will briefly discuss collider and 
gravitational wave signals for our scenario in Section |8} after considering the issue of complex 
phases in regards to EWBG and EDM constraints in Section [7) 



6 A Benchmark Example 

In this section we perform a numerical analysis of the EWPT and DM properties, without some 
of the approximations used in the previous sections. For instance, we fully take into account 
small i?- violating terms {B/j,, Mi and M2, though not A-terms), and various small vevs. We also 
include the effects of the temperature-dependent singlet terms [i.e. those proportional to C52 and 



cs in Eq. (27)]. However, we do not perform a full effective potential analysis, restricting rather 



to the thermal mass contributions at leading order. We also do not include CP-violating phases 
(see the next section for further discussion on this point). 

For a given choice of parameters, we flrst minimize numerically the potential at zero tem- 
perature, and determine the spectrum of Higgses, neutralinos and charginos, as well as their 
composition. Only the dimensionless ratios of the dimensionful input parameters are meaning- 
ful, since we can always rescale them a posteriori to normalize to f = 174 GeV. We then consider 
the finite-temperature potential to find the critical temperature. In order to do so, we increase 
T, and look for all real extrema of the potential (which is just polynomial in the approximation 
we are using), identifying EW-preserving (where only S has a vev) and violating minima. We 
then compare the potential energies associated with these minima in order to identify the global 
one, for the given T. By iteration, we can then identify the critical temperature, deflned by 
degeneracy of the EW symmetry preserving and violating minima. Note also that since we can 
take all squark and slepton soft masses to be positive and large (unlike in the MSSM, they do 
not play a role whatsoever in making the transition flrst-order), we do not need to worry about 
color- or charge-violating minima. 

As an example, consider the set of parameters (in GeV units where relevant) given by 
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We also take g ~ 0.65, g' ~ 0.35 and yt ~ 1, which results in c^ ~ 0.92, C52 ^ 0.24 and 
Cs ~ —2.78 GeV, from Eqs. (26) and (27). For the zero-temperature radiative corrections 
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we use MsusY = 2 TeV. The zero-temperature Higgs vevs are given by Vu ~ 173.99 GeV, 
Vs ^ 62 GeV, Vd ~ 1.2 GeV, w^ ~ -7 x 10"^ GeV. Note that Vu/vd ~ 150, illustrating that 
in this framework "tan/3" can be much larger than in the MSSM context (see next section). 
However, this ratio can easily be changed by changing 6, which barely has an effect on the main 
physical properties of the model. Also, the triplet vev is sufficiently small to be irrelevant from 
the point of view of EW precision constraints. 

The spectrum of CP-even (mj^J, CP-odd (^^aJ and charged {m^±) Higgses, in GeV, is 
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while the neutralino and chargino spectra are given by 
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It also of interest to note the composition of the two lightest neutral CP-even Higgses: 



Hi 

H2 



0.88 hi - 0.003 h^d + 0.48 s - 0.003 Tl 



0.47 hi - 0.008 hi, - 0.88 s + 0.005 Tf 



R ! 



(52) 
(53) 



and of the LSP: 



Xi 



° -' 0.676 



0. 12 w^ + 0.05 m + 0.35 T^ - 0.54 ^ - 0.35 H" 



(54) 



The content of X2 is very similar due to the pseudo-Dirac nature of the neutralinos. For this 
parameter point, we find that the relic density is Q^oh"^ ^ 0.11, in accord with the WMAP 
constraint |l5], while the LSP spin- independent cross-section is given by o-^o^_>^o^ ^ 4.5 x 
10"^^ cm^, somewhat below the current XENONIOO limit [i3] . 

Turning now to the finite-temperature analysis, we find that the critical temperature is 
Tc ^ 71.0 GeV, while Vc/T^ ~ 1.34, exhibiting a strongly first-order EWPT, even when taking 
into account the points made in [57J. For comparison, the analytic formulas given in Section 111 
give T^anaiytic ^ g2.5 GeV, while (wc/Tc)''"^^^*^'^ ^ 1.86. The difference arises from the temperature- 
dependent singlet terms proportional to 05-2 and cs- Nevertheless, we see that the physics is 
correctly captured by the simplified analysis. These results will shift slightly under a more 
detailed full effective potential analysis, but we do not expect that the conclusions will radically 
change. 
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7 QP Phases - EWBG and EDM's 

In Sections |3] and |4] we studied the scalar potential in the case that all the relevant parameters are 
real. This allows for a transparent understanding of the mechanism behind a strongly first-order 
phase transition, and its connection to the DM sector, that is typical in the scenarios we consider, 
as discussed in the previous sections. However, the production of a baryon asymmetry during 
the EWPT requires CP- violating phases. A complete treatment of this question, namely a study 
of the CP- violating sources that enter in the transport equations, their diffusion in front of the 
bubble wall, the communication of CP-violation to the LH quark sector and the subsequent 
processing by sphalerons, is beyond the scope of this work. Nevertheless, we will make a few 
remarks in connection to the production of the BAU that suggest that a large enough baryon 
asymmetry can be induced, while being easily consistent with the currently null EDM searches. 
To be definite, we will frame our discussion in the context of "approach IV" described in 
Section [21 and discussed in more detail in Appendix [B] We consider initially the exact \J{\)r 
symmetric limit. The model contains four Higgs doublets {11^-, Hd-, H'^, H'^, of which only H^ 
and H'^ acquire non-zero vevs), a SM singlet S, a SU(2)l triplet T, and a SU{?>)c octet O. The 



superpotential and the soft breaking terms are given by Eqs. (64) and (65) of Appendix Kl to 
which one adds the "supersoft" operators [see Eq. ([5|] iMoa^a^a + V^MoaDa^a + h.c, where 
Ai = B, \2 = W, A3 = g, T,i = S, T,2 = T and S3 = O, and the Da are the auxiliary fields 
in the SM vector superfields. For simplicity, we have focused on a limit where H'^ and H'^ are 
somewhat heavier than the weak scale (due to a large /i'-term) and can effectively be integrated 
outj^ However, we will keep them below initially in order to count all possible phases more 
transparently. 

The physical phases in the model, associated to the Higgs sector, are given by 

Arg(t,M^J , ATg{B,M*^\) , Aig (BtM*^^) , Arg (SoM*^) , (55) 

Arg(A5/i^) , Arg(AT/i^) , Arg (Al/i'^6') , Arg (A^/i'^6') . (56) 

It is convenient to choose M^^ to be real (by redefining the S, T and O superfields). Also, the 
phase of b' can be absorbed in H'^, for instance, and then we can choose /i and fi' to be real by 
rephasing the Hd and H'^ superfields. Thus, we can choose the physical phases to reside in A^, A'^, 
At, A'j., ts, Bs, Bt and Bq- The octet fields play no role during the phase transition, so we can 

^^Such a limit need not be essential to reach our conclusions: for instance, even if the primed Higgses have 
masses at the EW scale, as long as their vevs are somewhat suppressed our analysis of the phase transition 
should remain as a reasonable approximation. As {H'^ becomes larger (note that {H'^) can remain small due 
to the approximate i?-symmetry), the analysis will become more complicated, but we don't think the fact that 
the system can easily display a strongly first-order EWPT will change. Similarly, a connection to DM should 
remain, as well as the suppression of EDM's, as discussed below. 
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set Bq aside. Also, since the triplet vev is required to be small by EW precision constraints, any 
possible phase in Bt is likely to play only a minor role. In addition, given that {Ha) = {H'^) = 
(i.e. we take their squared soft masses to be sufficiently positive), and that we can integrate out 
the heavy H'^ and H'^, the Higgs potential can be seen to depend only on HlHu, H^Hd and S 
(setting T ;^ 0, together with the vanishing of the octet O and squark/slepton fields). It follows 
that non- vanishing phases in the microscopic parameters can only induce a phase in the singlet 
vev, while all the Higgs doublet vevs either vanish, or can be taken to be real by an SU{2)l 
transformation. 



The phases of A^ and A^ enter in the chargino and neutralino mass matrices [see Eqs. (21) 



and (22)], as does a possible phase in {S) (A^ and A'^. do not enter, as long as the H'^ and 
H'^ superfields are somewhat heavy). The S vev gives an interesting effect, not present in the 
MSSM, when Arg(S') is spacetime-dependent. Such a case was considered in Ref. [H] in the 
context of the related model of Ref. [19] (which does not have the U{1)r symmetry, but whose 
tree-level potential has the same form as ours). For the purpose of generating the BAU, the 
most important source of CP-violation comes from the chargino sector, at second order in the 
gradient expansion (assuming no degeneracy). In Ref. jH], it was found that the effect of the 
spacetime-dependent singlet complex vev easily leads to a significant baryon-to-entropy ratio 
V = ns/s oc A9s/{lwTc), where /^ is the bubble wall thickness (typically large compared to 
1/T ~ ^/Tc), and A^^ is the change in the phase of the singlet across the bubble wall. Thus, 
the phases in tg and/or Bs can generate a significant baryon asymmetry within the framework 
via a chargino source with CP-violation arising from the change in the singlet phase, even if Ag 
and At are real. 

As emphasized in the Introduction, some level of R-violation is generically expected once the 
gravitino mass is generated. In fact, we pointed out in Section |5] that small Majorana gaugino 
masses can play a crucial role in order for the LSP to fully account for the DM content of the 
universe. This suggests a second, qualitatively different possibility for generating the BAU. If 
all the phases in the Higgs potential vanish (as assumed in the analysis of Section 111), while the 
CP-violating phases arise only from the suppressed Majorana gaugino masses. Ma, one expects 
that the CP-violating sources will be proportional to Im(Ma)/MD^. Such a suppression might be 
welcome in regions where a too large baryon asymmetry is produced, as suggested by the results 
of IH] . It would be interesting to further study the detailed aspects of EWBG in the previous 
scenarios (note that the model studied in [H] required sizable gaugino Majorana masses). 

The above mechanisms for generating CP-violating sources during the EWPT should remain 
operative in the context of "approach I" described in Section [2j However, there can exist 
differences between approaches I and IV when it comes to EDMs. In order to see this, it will be 
useful to ffist comment on "tan /3" . We have assumed that the EW vev, v, is carried mostly by one 
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Higgs doublet (-ff°), and therefore in a sense we are always considering a "large tan/3" scenario. 
As is well-known, in the MSSM one often finds observables that are tan/3 enhanced/suppressed. 
In particular, large tan/3 enhancements arise in the down sector when the vev that gives rise 
to the down-type fermion masses is much smaller than v. Such enhancements are, however, 
tied to how these masses are generated, as illustrated by the two operators in Eq. Q. If the 
down-fermion masses arise from the first operator (as in approach I with a small i?-violating 
Bfi-term being responsible for nid and me), the situation is "MSSM-like" in regards to such 
enhancements. However, if the down-fermion masses arise from the second operator (as in 
approach I with the Dobrescu-Fox mechanism [21], or as in approach IV), then the vev of H^ is 
not related to the down Yukawa matrices, yd,e- Instead, the measured fermion masses always 
relate the Yukawa couplings to Vu ~ v, and the source of tan /3 enhancements disappears. Notice 
that, as in approach IV, the down masses may be ultimately connected to a small vev like v'^ 
(see Appendix [B]). Thus, one can define two different (large) vev ratios, Vu/vd and Vu/v'^. One 
may then wonder if the second ratio can play the role of the MSSM tan/3, but as we will see the 
i?-symmetry prevents "tan /3 enhanced" terms from appearing. 

The above remarks will be useful when considering the issue of EDM's, which we will first 
discuss in the context of approach IV. As emphasized in [21], the approximate U{1)r symmetry 
leads to a significant relaxation of the constraints from the electron and neutron EDM's (see [15] 
for a review). The point is that, at one-loop order, the EDM's induced by a squark/slepton 
and gluino/chargino/neutralino loop require LR mass mixing in the sfermion sector. In the 
MSSM and many variants, including the model of Ref. [I9], these LR mixings arise from A- 
terms and the /x-term, the latter effect being tan/3 enhanced in the down-type sector. In our 
scenario the situation is different: the A-terms are at least loop-suppressed (if they arise from 
anomaly mediation), while the U{1)r symmetry forbids the "usual /i-term", i.e. a superpotential 
term coupling the Higgses responsible for the up and down-type fermion masses (in our case 
HuH'^. The "/i-terms" allowed by the U{1)r symmetry, ftHuHa and ^H'^H'^, do not contribute 
to LR sfermion mixing, as long as H^ and H'^ do not acquire vevs. This holds even in the 
presence of (small) Majorana gaugino masses, as required by the DM relic abundance (discussed 



in Subsection 5.1). Note that such Majorana gaugino masses will induce, at one-loop order, a 
b^HuHd term, with 6^ ~ [aaMafi/{27i)] log(A//i) [29j where a^ is a SM fine-structure constant 
and A is a UV cutoff. Such a 6-term induces in turn a very small /^-violating vev for {Hd) ~ 
Vvb^/f^'\j 1 that translates into very small LR sfermion mixings in the down sector. In addition, 
since the down-type fermions get their masses mostly from the H'^ vev, instead of Hd-, such LR 
mixings do not carry any "tan /3" enhancements. Thus, we conclude that there are no constraints 
from one-loop level EDMs in our scenario, even if the sfermion masses are light (~ 1 TeV). This 
can be contrasted to the case of pseudo-Dirac gauginos without the U{1)r symmetry, studied 
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in [HH]: while the Ma/Mj:)^ suppression can be common to the generation of the BAU and EDM's, 
in the f/(l)_R models the additional LR suppression renders the associated EDM's completely 



harmless, ^^ Similarly, two-loop Barr-Zee-type diagrams [59], involving a chargino in the loop, 
are also expected to give a contribution to EDM's suppressed by Ma/M^^, and are therefore 
negligible. As pointed out in [21], the leading contribution to EDM's in f/(l)/? symmetric 
scenarios corresponds to the Weinberg three-gluon operator [60], {w/3)G ■ G ■ G, which could 
be observable in the near future if the relevant CP-violating phases are order-one. However, 
these phases are associated with the gluino/octet sector, and therefore are not closely connected 
to the relevant phases for EWBG. Also, if the relevant CP-violating phases responsible for the 
BAU arise from the Higgs sector, as discussed in the first EWBG scenario above (say with very 
small/real Majorana gaugino masses), the CP- violation is communicated to the SM fermion 
sector only at a high-loop order, hence it is unconstrained by EDM bounds. 

We end this section by commenting on EDM's within "approach I". If the down-type fermion 
masses arise at one- loop order as in ^29j-^° then the situation is similar to approach IV described 
above, since the down-fermion masses effectively arise from the second operator in Eq. Q [see 
discussion on tan/3 above]. If, on the other hand, the down-type fermion masses arise from a 
small Hfi vev induced by a i?- violating S/i-term, then one can expect large tan /3 enhancements 
in the one-loop induced EDMs, that may compensate the Ma/M^,^ suppression associated with 
the pseudo-Dirac nature of the gauginos. A more detailed study is then necessary to estimate 
the bounds on the sfermion masses, but one can nevertheless expect an improvement compared 
to the situation in the MSSM (see e.g. [61]). This conclusion also holds for EDM's induced at 
2-loop order. 

8 Other Experimental Signatures 



As mentioned at the end of Subsection 5.2, a combined set of measurements from a variety 
of experiments in particle physics and astrophysics will be needed in order to fully test the 
framework and zoom in on the preferred region of parameter space. In addition to DM direct and 
indirect detection, or EDM measurements, this would include collider physics signals at the LHC 
and other future collider experiments, as well as a possible gravity wave signal arising from the 



^^We note that the pseudo-Dhac nature of the gaugmos by itself may be sufficient to aUow for a successful 
EWBG, while being consistent with EDM searches, for natural values of the sfermion masses: the singlet can 
be responsible for a strongly first-order phase transition, such that in spite of the Ma/M]j^ suppression of the 
CP-violating sources, a sizable BAU can be generated. At the same time, the 1-loop EDM's are suppressed by 
the same Ma/Mn^ factor, allowing for lighter sfermions than in the MSSM. 

^°However, note that due to the suppression in the Majorana gaugino masses, which comes on top of the one- 
loop suppression, it may be necessary to take rather large values for /i in order to generate the bottom Yukawa 
coupling 29 . Thus, this scenario may require some fine-tuning in the Higgs sector. 
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strongly first-order EWPT. Although a comprehensive study of phenomenological consequences 
will have some model-dependence and is beyond the scope of this paper, we will attempt to 
outline the characteristic signatures of the framework which depend only on its crucial features. 

8.1 Collider Signatures 

The framework has a number of interesting collider signatures which can be probed at the LHC 
Some of the signatures arise as a consequence of the broad features of the framework considered 
in this paper, while others arise as general consequences of the R symmetry and hence share 
signatures with other i?-symmetric models considered in the literature. We will focus primarily 
on the former and only briefiy mention the latter, directing the reader to the relevant papers for 
reference. 

The class of models considered here have a characteristic spectrum. The lightest degrees of 
freedom consist of the lightest CP-even Higgs, and the lightest chargino and neutralino, with 
the lightest neutralino being the LSP. They all have masses < 120 GeV. The other Higgses, 
charginos and neutralinos generically are heavier, spanning a large range between around 150 
to several hundred GeV. However, in order to satisfy the constraints on the T-parameter, the 
triplet scalars have to be quite heavy, in the multi-TeV range. What about squarks, sleptons 
and gluinos? Since the physics underlying electroweak baryogenesis and DM in this work is 
essentially independent of them, and given the lack of significant constraints from EDM bounds, 
the masses of these particles are quite unconstrained, and can range from a few hundred GeV 
to multi-TeV. 

The lightest CP-even Higgs in the allowed region of parameter space generically has a non- 
trivial singlet component. Therefore, the LEP constraints on the Higgs mass are somewhat 
relaxed, the exact amount depending on the parameters in detail. Although there exist additional 
tree-level contributions to the Higgs quartic in a model with the coupling X^SHuHd, the extra 
contribution vanishes in the f ^ — ?■ limit, implying that radiative corrections are important in 
raising the Higgs mass. It is more challenging to discover the lightest CP-even Higgs in this 
class of models than in the SM or in the MSSM, due to two main reasons. First, as mentioned 
above, the Higgs has a sizeable singlet component which reduces its coupling to gauge bosons 
and quarks. Second, the decay mode h ^ XiXi is generically available, allowing the Higgs to 
decay invisibly with an appreciable branching ratio, ^^ In such a case, vector boson fusion is 
expected to provide the most effective search channel at the LHC j62|. On the other hand, the 
orthogonal, and heavier combination of h^ and s (call it H), can have a sizable coupling to gauge 
bosons and to the top quark, and be in a range where its decay into W/Z pairs has a sizable 



"^The lightest Higgs is expected to be below the weak boson production threshold. 
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Figure 6: Correlation between the X2 ^^^^ frame decay length and the xl relic abundance in the model of 
Section [6J as the U{1)y Majorana gaugino mass, Mi, is varied in the range [1, 10] GeV. The dashed horizontal 
line indicates approximately the minimum decay length that is measurable. 



branching fraction. Such a state may be looked for m gg ^^ H ^ WW/ZZ. It hardly needs 
to be emphasized that Higgs physics in this framework is extremely rich and interesting, and 
should be thoroughly investigated. 

The signatures related to the lightest neutralino and chargino are also very important since 
they are intricately tied to the DM physics. The lightest chargino in the allowed region tends 
to be close to the LEP bound; hence should be discovered at the LHC. What about the LSP? 
The pseudo-Dirac nature of the LSP gives rise to two quasi-degenerate states x? and X2 ^plit by 
a small amount (5m, causing X2 ^^ decay to Xi- Since 5m <^ m^o, the decay length L could be 
large in some cases, leading to a displaced vertex. However, since L ~ 6m~^, the decay length 
is quite sensitive to 6m and could span a large range. 

Note that X2 decays to Xi and SM fermion pairs (//) through an operator which is precisely 
the same as the one relevant for computing the relic abundance of the LSP, viz. Xi and Xo, co- 
annihilation to // through Z-exchange! As observed in [63], if the X2 decay length is measurable, 
it allows for a non-trivial "measurement" of the relic density, solely based on collider data, that 
can then be compared to the cosmological observations. The point is that the X2 decay width, 
given by 



nx'2 -> X?//) 



C"4 Sm^ 
120^ M| 



(57) 



,0 and on a "coupling constant" C which can be effectively 



depends only on 6m, = m^o — m. 

identified with the constant C appearing in the annihilation cross-section, Eq. (46). The only 

difference is that the bottom quark contributes to C, but not to C", which results in C ~ 1.18C". 
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One can then see that there is a correlation between Q^oh"^ and the X2 decay length (given the 
two measurable quantities m^o and 6m) that is independent of the detailed composition of the 
LSP. A difference between our situation and that of Ref. [63] is that in our case the annihilation 
proceeds in the vicinity of the Z resonance, while in [63] both the annihilation and the decay 
width of the heavier pseudo-Dirac state were assumed to proceed via a very heavy state, leading 
to an effective contact interaction. As a result, we find that it is possible for the decay length 
to be macroscopic for pseudo-Dirac mass splittings that result in an order one relic density. 

As an example, we show in Fig. M the X2 decay length versus Vt^ah'^ in the benchmark 
example of Section [6] We only vary the Majorana gaugino masses, imposing for concreteness 
the "unified-like" relation M2 = 2Mi, and letting 1 GeV < Mi < 10 GeV. In this range, the 
mass sphtting between the two lightest neutralino states (with masses of about 60 GeV) varies as 
0.85 GeV < 5m < 10 GeV. We then find that the X2 decay length in the rest frame is in the range 
10 cm > L > 0.5 yum, and the rehc density satisfies 0.05 < Vt^ph? < 0.19. Other properties, such 
as the strength of the phase transition, are not sensitive to the Majorana masses. Also, in the 
full range above, the spin-independent cross-section is consistent with current bounds. Thus, we 
see that 0{1) relic abundances can be associated with observable (> 10~^ mm [HUES]) displaced 
vertices. Note that if the LSP is only one of several components of the total DM abundance, 
the above method (provided L is measurable) can give a very useful handle in determining its 
abundance. However, it is also possible that X2 is stable on collider time scales, in which case it 
could not be distinguished from the real LSP in a collider environment. 

Before moving on to other signatures, we briefly mention some collider signatures which 
depend only on the existence of an approximate i?-symmetry and are shared with other R- 
symmetric models in the literature. In the pure i?-symmetric limit, it can be shown that 
cross-sections for many scattering processes, such as qq' — >■ {(JlQ'l: QrQ'r} : I'i' ~^ {(lL(j'R,(JR(j'L}y 



Q9 ^ {qLgD,qRgD}, and gg -^ [gogD^gBgD]^ vanisfQ p)]. Thus, these scattering processes 
are suppressed if the /2-breaking is small. The suppression of Majorana masses also implies 
suppression of same-sign lepton signals. In the Higgs sector, the presence of the i?-charge 2 
scalar Hd implies that it can only be pair-produced at colliders. At the LHC, Drell-Yan pro- 
duction mediated by weak gauge bosons is likely to be the dominant channel [67]. Also, in 
the i?-symmetric limit each Hd decays to two neutralinos (with i?-charge 1) giving rise to four 
neutralinos for events with pair-produced HdS. Finally, since the i?-symmetric SSM can lead to 
a natural suppression of FCNCs even for 0{1) flavor violating soft masses, the flavor violations 
in the soft mass-matrices can show up in variety of ways at the LHC [2l], which is worthy of 
significant study. 



^cjD stands for the Dirac gluino. 
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8.2 Gravitational Waves 

A strongly first-order EWPT proceeds by formation of bubbles of the broken phase which expand 
into the unbroken phase. Near the end of the phase transition, bubbles of the broken phase collide 
with each other, breaking the spherical symmetry and leading to production of gravitational 
waves [68]- [72]. Thus, a strongly first-order EWPT can give rise to gravitational waves, which 
are in principle observable by space-based interferometers like LISA or BBO. 

A detailed study of the feasibility of observing a gravitational wave signal arising from the 
EWPT within our framework is beyond the scope of this paper, and we will limit ourselves to a 
few remarks. As explained in Section |4| the structure of the finite temperature potential within 
the framework is qualitatively different than in the MSSM or the SM, which naturally allows the 
possibility of a strongly first-order EWPT. The crucial feature is the presence of a barrier due 
to a negative effective quartic term in the effective potential, which is balanced by a "tower of 



operators" arising at higher order (see Eqs. (31 )-(33) and the discussion below). It turns out that 



a stronger phase transition proceeds at a lower temperature and creates larger bubbles, shifting 
the peak of the gravity wave spectrum to lower frequencies, below the best sensitivity range 
of satellite experiments like LISA and BBO. However, as explained in [TD], since the gravity 
wave spectrum has a much milder fall-off with frequency (~ /~^°) than thought earlier [73], it 
is expected that BBO would be sensitive to gravity waves from the EWPT. Also, as pointed 
out recently in [71], EWBG and gravity waves depend on different velocities associated with 
the expanding bubble wall. While gravity waves depend on the wall velocity Ku as a whole, 
electroweak baryogenesis is sensitive to the relative velocity f + between the wall and the plasma 
in front, which is in general smaller than Vw, the difference being more pronounced for stronger 
phase transitions. Thus, it may be possible to have an observable gravity wave signal from BBO 
or LISA from a strong electroweak phase transition responsible for electroweak baryogenesis j74j . 

9 Summary and Conclusions 

In this work we have considered a supersymmetric extension of the SM that have an approximate 
U{l)fi symmetry. We point out that this class of models can naturally lead to a strongly first- 
order electroweak phase transition, and possibly to the successful production of the observed 
baryon asymmetry, without the tensions that plague the MSSM. We have further pointed out 
that there is a close connection between the EWPT and the properties of DM (see [75] for a 
study of the DM/EWBG connection in the context of the MSSM). 

The basic observations derive from the (pseudo) Dirac nature of gauginos in i?-symmetric 
models. This is motivated by the fact that Majorana gaugino masses are forbidden by the R- 
symmetry, which leads one to introduce superfields in the adjoint representation of the SM gauge 
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group. In particular, the Dirac partner of the bino arises from a SM singlet superfield, whose 
scalar component plays an essential role in leading to a first-order EWPT. Since this is a tree- 
level effect, the phase transition is easily strong, even for Higgs masses comfortably above the 
Higgs LEP bound, thus creating the necessary conditions for generating the BAU. We find that 
the interesting region of parameter space is also characterized by a sizable Higgsino component 
of the lightest neutralino (the LSP), which results in an appealing DM/EWBG connection. 
This fact, together with the (pseudo) Dirac nature of the LSP, leads to rather interesting DM 
physics: at the time of freeze-out, the LSP behaves similarly to a Dirac fermion, and has a 
sizable annihilation cross-section into fermion pairs via Z-exchange. This cross-section can be 
modulated by a small Majorana mass splitting. On the other hand, for annihilations today, or 
in scattering against nuclei, the LSP behaves like a Majorana fermion, thus being consistent 
with current constraints. Nevertheless, one expects that such a DM candidate will be observable 
in upcoming direct detection experiments, and possibly also at neutrino telescopes, such as 
IceCube/DeepCore. 

We have further stressed that, unlike in the MSSM, the CP-violation that is relevant in the 
generation of the BAU is typically not constrained by EDM searches (although EDM signals 
are possible within the framework). We have also commented on the collider prospects, which 
are characterized by a rich Higgs sector, and a possibly very interesting signal associated with 
the two semi-degenerate lightest neutralino states - in the form of a visible displaced vertex 
in association with missing energy, arising from the decay of X2 ^'^ Xi- Moreover, the collider 
measurements of the decay length L, m^o and 6m allows the "measurement" of a cosmological 
observable - the relic abundance of the LSP. Due to the strength of the EWPT, gravitational 
wave signals represent another exciting possibility. 

In conclusion, the present class of models represents a well-motivated extension of the SM 
(addressing the hierarchy problem, while significantly alleviating the SUSY fiavor and CP prob- 
lems), that can also easily lead to a successful generation of the BAU during the EWPT, while 
providing a non-standard DM candidate. In all these regards, it compares favorably to the 
MSSM. More detailed studies of the points made above are worth pursuing and are left for 
future work. 
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A Expressions including i?-breaking and vt ^^ 



Here we collect the relevant expressions when i?-symnietry breaking terms by Majorana gaugino 
masses and the h = Bfi term are included, and when vt y^ 0. However, we still assume that 
74-terms are negligible. Now, unlike in the analysis in Section [3} the H^ degrees of freedom no 
longer decouple from the system and a small vev v^ is induced for H^. Also, the triplet vev vt will 
be small but non-zero. For simplicity, we will only consider a vacuum which is CP-preserving, 
i.e. only the real parts of the fields get vevs. The minimization conditions are then given by: 
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We now write the components of the general four-dimensional neutral CP-even, neutral CP- 
odd and the charged Higgs scalar mass-squared matrices. Since these are symmetric, there are 
ten independent components. For the CP-even scalars in the (/i°, /i°, Sr, T^) basis, and denoting 
the mass-squared matrix as A^^, one has: 
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For the CP-odd scalars in the [G^,A^,Sj,Tf) basis, and denoting the mass-squared matrix as 
Ai\, one has 
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The mass-squared matrix elements for the charged Higgs scalars, written as l-L^ Ai^^^l-L , with 
n+ = {H^,HJ,T^,T^} and H- = {R- ,H, ,T- ,T^}, are: 



41 



M' 



H±A1 



9_ 

4 



-J2 



-C2/3 + si AA 



+ {XsVs - XtVtY + V2 (g'MD.Vs + qMd^vt) + m]j^ , 



M' 



M' 



H±,21 

2 
H±,31 



b+ ( |--A^ + A|)s/3C^t;2 , 



S/3 



A^2 



//±,41 



Al^ 



_ff±,22 



M' 



W 



H±,32 

2 
H±A2 



2gMD^ + V2g^VT - 2V2\t{\sVs - \tVt) 
2gMD, + V2g^VT + 2V2 \t{\sVs + XtVt) 

v^ + {XsVs - XtVt^ - V2{g'MD,v, + gMo^vr) + m]j^ , 



S/3 






J 2 



A^2 






H±,33 

2 
fl"±,43 

2 
H±AA 



2gMD, + 72 ^V + 2^2 AT(A,t;, + Ayt^r) 
- -2gMD2 + V2g^VT-2V2\Ti\sVs-\TVT) 

( |-C2/3 + 2A| sUv^ + g^ vl + 4^2 ^ Md^vt + m^ + 2 M 

-(72 4 + 25t + 2M|,^ , 

(-|-C2/3 + 2A| 4 j t;2 + g^vl- AV2gMD,VT + m? + 2 M 



2 
D2 ' 



2 

1)2 



Similarly, for charginos in the {T+, PF"*", i/+ ; T , PF , -^^ } basis, one has: 
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Finally, the neutralinos become Majorana in character once /^-breaking is included. Thus, in 
the {i-B, iW^, H^, T, S, H^} basis, the neutralino mass matrix is given by: 
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B Approach IV: Details of the Model 

In this section, we explain in some detail approach IV discussed in Section |2j In particular, in 
addition to the fields in the _R-symmetric SSM, we include a vector-like pair of SU{2)l doublets 
H'^ and H'^ with R[H'^] = 2,R[H'^] = and Y[H'J = l/2,Y[H'a] = -1/2. The i?-charges of all 
the quark and lepton superfields are equal to unity, while R[Hu] = and R[Hd] = 2. This allows 
the following superpotential consistent with the gauge and i?-symmetries: 

{y'^QD^ H', + f^' HiH', + KS KH', + X'^T KH',) . (64) 

The soft supersymmetry breaking terms in the potential consistent with the U{1)r symmetry 
are given by: 

^soft = mljH^\^ + mj,jHd\^ + ml^\H'J^ + ml,jH'f + {b'HuK + h.c.) + 

{t,S + B,S^ + Bt T^ + h.c.) + ^ ml^ |Qi|2 + ^ m^^ \U\^ , (65) 

i i 

where, for simplicity, we omit the octet scalar terms and a possible cubic term for the singlet 
(expected to be small). We denote all the squarks by Q and all the sleptons by L. We assume 
that all the soft mass parameters {m'jj^,m'jj^,m'jj, ,771"^, ,mQ,m'j^,b', Bs,Bt} are parametrically 
of C(Mg^Q£^) -close to the TeV scale. We will also assume that /i ~ Mgoft, but that /x' is larger 
than Msoft- Therefore, the scalar and fermionic components of the H'^ and H'^ superfields are 
heavy (with masses ~ /i') and decouple from the EW scale physics. 

Due to the assumption that /i' ^ Msoft, H'^ does not get a vev, and neither does H^ if 



mjj > 0. On the other hand, the b' term in Eq. (65) induces a small vev for H'^ once H^ gets a 



non-zero vev. This vev will then induce masses for down-type fermions: 

/ Vub' , Vub' 

The above result can also be seen in terms of operators by integrating out the heavy superfields 
H'^ and H'^. Since /z' is assumed to be larger than Msoft, it is possible to integrate out the entire 
supermultiplet. Solving the equation of motion for the superfield H'J gives: 

H',^ = --r^, (y'd ^"Q T^cD^ + y'e T^'^L V^E'^) + ■■■ (67) 

4/i' 

Substituting the above solution for H'^^ in the b' term, written as jdl^O [X^ X / Ml)HuH'^ + h.c. 
gives rise to: 
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and therefore 



, VuM^^ 



2 
soft 



^ ~ yi^^^br^- (6* 



Mf,. 



Note that due to the suppression in v'^ from //' ^ Mgoft, the bottom Yukawa couphng may 
be of order one, corresponding to a large tan /3 = f u/f ^ scenario. However, as a result of the 
structure of the f/(l)/j symmetric model, this does not necessarily lead to the typical large tan/3 
enhancements that characterize certain observables within the MSSM. An example occurs for 
the electron and neutron EDM's, as described in Section [7j 
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